A Three Dimensional Kasteleyn Transition: Spin Ice in a [100] Field 
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We examine the statistical mechanics of spin-ice materials with a [100] magnetic field. We show 
that the approach to saturated magnetisation is, in the low-temperature limit, an example of a 
3d Kasteleyn transition, which is topological in the sense that magnetisation is changed only by 
excitations that span the entire system. We study the transition analytically and using a Monte 
Carlo cluster algorithm, and compare our results with recent data from experiments on Dy2Ti207. 

PACS numbers: 75.10.Hk, 75.40.Cx, 75.40.Mg 



Topological phases occur in connection with new liquid 
states and concurrent exotic transitions out of them. A 
central idea is the emergence of degrees of freedom which 
cannot be defined or manipulated locally. Their presence 
has profound ramifications for the static and dynamic 
properties of such phases. They can exhibit algebraic 
correlations in the absence of criticality [ij] and unusually 
slow relaxation towards thermal equilibrium [3, Q . 

A magnetic system whose topological properties have 
a particularly simple and transparent origin is known as 
spin ice [1, Q (Fig. [T]) : in its macroscopically degenerate 
ground state, a magnetic version of the ice rules apply. 
These stipulate that for each tetrahedron of the corner- 
sharing pyrochlore lattice, two of the four spins point in, 
and two point out of the tetrahedron along the local body 
centred diagonal axes. As a result, the lattice divergence 
of the spin field vanishes everywhere: V • = 0. As a 
consequence of this local constraint, the magnetisation 
of each (100) plane is the same Q and is therefore a 
topological quantity which can be changed only by mak- 
ing a simultaneous change throughout the system: the 
smallest such magnetic excitation involves a set of spins 
on a string spanning the system. An attractive feature 
of spin-ice materials is that one can couple directly to 
this topological quantity using a uniform magnetic field, 
which lifts the macroscopic degeneracy of the zero-field 
ground states. 

In this paper, we focus on the behaviour near satura- 
tion, when the ordered g = state [1, 0] is reached by 
applying a magnetic field along the [100] direction as il- 
lustrated in Fig. [TJ In the saturated state at large fields, 
the ice rules prohibit local fiuctuations in the limit of low 
temperature T. The resulting phase transition as the 
field is lowered is a three dimensional (3d) example of 
a Kasteleyn transition Q, first proposed in the context 
of 2d dimer models Its hallmark is an asymmetric 
character, appearing to be first-order on one side, and 
continuous on the other. The Kasteleyn transition has 
in the past been associated with soft matter: the trans- 
gauche disordering transition for polymers embedded in 




FIG. 1: (Color online) Left: Pyrochlore lattice showing g = 
spin structure (black narrow spins) and a string defect (red 
thick spins). Right: Bethe lattice with a central tetrahedron 
(red) and first layer of nearest-neighbour tetrahedra (blue). 



lipid bilayer system, observed via density measurements 
in dipalmitoyl lecithin (DPL) [l3|. To our knowledge 
this is the first example of a three-dimensional Kasteleyn 
transition in a magnetic system. 

We analyse the transition using a combination of ana- 
lytical and numerical tools. By mapping the strings onto 
bosonic worldlines, we obtain an analytical expression 
for thermodynamic quantities close to the critical point, 
and for correlation functions at all values of the magnetic 
field. These predictions are confirmed to high precision 
by numerical simulation using a Monte Carlo algorithm 
incorporating a non-local cluster update. In addition, 
we find that a calculation on a Bethe lattice provides an 
excellent description of the thermodynamics of the mag- 
net away from the critical point. A crucial characteristic 
of this transition is that it is symmetry-sustaining, as 
recognised in Ref. 



11| : the first-order liquid-gas transi- 



tion reported there, however, is an artefact of the local 
dynamics employed at that time. We compare our results 
with magnetisation measurements on the spin ice com- 
pound Dy2Ti2 07 [13]: they represent well the observed 
behaviour at large field and temperature, but for lower 
values the experimental system falls out of equilibrium. 
We consider a spin ice model with nearest-neighbour 
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magnetic interactions only. The Hamiltonian reads 

(io) « 
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where Si is a vector of unit length, encoding an Ising 
spin constrained to lie along the local body centred cubic 
axes of the tetrahedron; J > is the ferromagnetic ex- 
change constant and h is an external field of strength h 
in reduced units, aligned along the [100] direction. This 
model is a simplification in that the real Hamiltonian 
is known to incorporate dipolar interactions of at least 
equal importance 13|, However, it has been demon- 
strated that a nearest-neighbour description is adequate 
as long as the ice rules are not violated, since corrections 
fall off with separation r as r^^ 15;]. Our results should 
therefore be robust for T <ti J, when behaviour depends 
only on the ratio h/T, whereas interactions between fur- 
ther neighbours may become important at higher T. 

Useful insights into the Kasteleyn transition come from 
discussing in detail the string excitations introduced 
above: start from the q = state, flip spins in succes- 
sive [100] layers to generate another configuration, and 
link flipped spins in adjacent layers with segments of the 
strings. This mapping between spin and string configura- 
tions is one-to-one provided we adopt a fixed convention 
for the paths followed by two strings through the same 
tetrahedron. Strings of finite length do not occur in ice- 
rule states as their ends cost exchange energy 0{J). In 
contrast, strings that span the system cost only Zeeman 
energy, of 2h/^/3 per segment. An isolated string also 
has an entropy, of In 2 per segment, and a finite density 
of strings can account for the macroscopic ground state 
entropy of spin ice. The total free energy for a single 
spanning string of L segments is 



G = L 



(2h/V3 
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From this one obtains the exact value of the Kaste- 
leyn transition temperature, Tk = 2/i/-\/31n2, and for 
h,T J, properties are a function only of T/Tk- For 
T < Tk the equilibrium phase is the q ~ state with- 
out any thermal fluctuations, while for T > Tk strings 
are present. Depending on their interactions, different 
scenarios are possible for the variation of their density, 
and hence sample magnetisation, with T/Tk- If strings 
attract, their density has a jump at the transition. At 
finite string density, if they form a commensurate crys- 
tal, there is a magnetisation plateau; and if they form a 
liquid, their magnetisation varies smoothly. In fact, by 
construction, strings have hard core repulsion, and the 
simulations we report below show only a liquid phase. 

By thinking of the strings as world lines for hard core 
bosons moving in two dimensions at zero temperature, 
with the [100] direction as imaginary time, we can ob- 
tain thermodynamic behaviour near the critical point. 



Under this mapping, the reduced temperature difference 
t = (T/Tk — 1) is proportional to the chemical poten- 
tial for bosons, the deviation of the magnetisation M per 
site from its saturation value Mgat is proportional to bo- 
son density, and the spin-ice free energy is proportional 
to the boson ground-state energy. The Kasteleyn transi- 
tion is between a vacuum state and a Bose condensate as 
the chemical potential changes sign. Here, three (2-1-1) 
is the upper critical dimension, since in the quantum de- 
scription the dynamical critical exponent takes the value 
z=2, a reflection of the fact that the [100] fleld renders 
the classical model highly anisotropic. Transcription of 
established results for the boson problem f\B\ gives for 
t > 



M = M,at[l-at-6tln(l/t)], 



(3) 



where a and b are constants. It follows from this that 
the heat capacity Ch and the differential susceptibility 
X diverge logarithmically in t as the critical point is ap- 
proached from the high-temperature side. Similar con- 
clusions for a dimer model were reached using ar gum ents 
based on the interaction between two strings [sI-IitI]. 

Correlation functions are of great interest throughout 
the partially magnetised phase, and at /i = are known 
to have an algebraic dipolar form that is a consequence 
of the ice rules [13] . The bosonic description of the sys- 
tem that we have introduced provides both a new way 
of viewing this result and also a means of finding how it 
changes with h. For a complete treatment, it is neces- 
sary to take into account the distinction between the four 
sites in the basis of the pyrochlore lattice when construct- 
ing the mapping between spins and strings. Moreover, a 
precise representation of the transfer matrix for strings 
in terms of a boson Hamiltonian would include multi- 
particle interactions, which are not expected to infiuence 
long-distance behaviour or critical properties, since both 
of these should be universal. We defer these details to a 
future publication and present here a simplified calcula- 
tion, in which bosons move on a square lattice and the 
imaginary time direction is approximated as continuous. 
We start from a boson Hamiltonian written in terms of 
spin S" = 1/2 operators as 



(4) 



where is an effective coupling constant and B oc h. In 
this incarnation, the Kasteleyn transition is a quantum 
phase transition in which the two-dimensional XY ferro- 
magnet is fully polarized by a field B coupled to the z- 
component of magnetisation, at a critical field strength 
B = Be- Moreover, M for spin ice is proportional to 
(iSf ) for the quantum magnet. We use a large S approx- 
imation to evaluate the connected correlation function, 
defined for spin ice as 



(5) 
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FIG. 2: (Color online) M vs. T /Tk obtained from simulations 
for the pyrochlore lattice (dots) and analytically on the Bethe 
lattice (solid lines) for h/J = 10"^(«), 0.13(a) and 0.58(B). 



FIG. 3: Ch and x vs t = (T/Tk — 1) obtained from simula- 
tions (•) and analytically (solid lines) for h = 10~'^J. The 
dashed line is a fit to a logarithmic singularity derived from 
Eq. ©. 



where the separation between lattice sites i and j is r 
within the (100) plane and z along the [100] direction. 
Since we have treated the lattice structure in a simplified 
way, our results apply only to intra-sublattice correla- 
tions. Writing B/Bc — cos 9, we find for large r and z 

, S- sin 6* 3 [z sin 61] 2 - (r 2 + [z sin 61] 2) 



(r2 + [zsin6']2)5/2 



For B : 
of Ref. 



the form of this expression reduces to that 
181 ] and has the cubic symmetry of the zero 
field problem. For finite field this symmetry is bro- 
ken and distances in the z direction acquire a scale fac- 
tor l/-\/l — {B/BcY- Close to the Kasteleyn transition, 
there is a second regime of behaviour, which also follows 
from Eq. Q but is not displayed in Eq. ([6|). It arises 
because correlations at distances shorter than the string 
separation are controlled by the behaviour of an isolated 
string, and hence reduce to the correlations along a 2d 
random walk, giving for example, C(0, z) ~ z~^. 

We have investigated the transition via Monte Carlo 
(MC) simulations, using a cluster algorithm 1^ 2^ 21). 
This allows us to avoid the apparent first-order disconti- 
nuities [m that arise from loss of ergodicity with single 
spin flip dynamics for T <C J. Each MC move consists 
of the reversal of all spins on a string, which in the limit 
T / J — > and with periodic boundary conditions must 
close on itself. The algorithm is efficient because in this 
limit all attempted flips are accepted below saturation. 
A typical simulation involves 4 x 10^ spins, studied for 
10^ MC steps at each temperature. Monte Carlo data 
for magnetisation (Fig. [2), differential susceptibility and 
heat capacity (Fig. [3]) show conclusively a 3(i Kasteleyn 
transition in the limit h <^ J . For T < Tk, the system is 
completely frozen, with M = Mgat and Ch = X = 0, while 
for T > Tk, M varies continuously with T/Tk- The be- 
haviour of Ch and x close to Tk, related by Ch — h'^x/T, 
is consistent with the logarithmic singularity expected 
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FIG. 4: (Color online) Finite-size scaling of spin correlations 
for hjhc = 0.33(B), 0.66(a) and 0.93(«), on log-log scale. 
Lines have slopes 3 and 1. 



from Eq. ([3]). Away from the limit h <^ J , there is no 
sharp transition: instead M varies smoothly with T . 

Data for the correlation function, illustrated in Fig. 01 
show the behaviour expected from Eq. ([5]). At the max- 
imum separation along the [100] direction in a system of 
linear size L (r=0 and z=L/2) correlations decrease with 
increasing L as , and increase in amplitude with h/hc- 
Close to the transition, a slower decrease is observed at 
short distances, consistent with the law expected for 
isolated strings. 

We supplement our simulations with analytical calcu- 
lations for a related model: a Bethe lattice (BL) of tctra- 
hedra, whose central element is illustrated in Fig. [TJ As 
on the pyrochlore lattice, each tetrahedron is connected 
to four others. Crucially, however, there are no closed 
loops of tetrahedra the statistical mechanics can 

then be solved exactly and a Kasteleyn transition ap- 
pears at Tk ■ Analytical results for the BL are very close 
to the simulation data for the pyrochlore lattice - in the 
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FIG. 5: Magnetisation M vs the internal magnetic strength 
fJ^oHint obtained analytically for the BL (solid lines) and ex- 
perimentally (points) for the compound Dy2Ti207 at Texp = 
1.8K; data from Ref. [11. 



absence of adjustable parameters (Fig.[2]and[3]) - provid- 
ing a posteriori justification for the BL approach. SmaU 
differences are apparent only for h <^ J: these arise be- 
cause the transition on the BL is necessarily mean field, 
and so does not have the logarithmic contributions to 
critical behaviour of Eq. ([3|) . 

What about the experimental situation? Studies of 
magnetization as a function of [100] field strength at 
fixed low temperature have been reported for the spin 
ice materials DyaTiaOr [l3| and Ho2Ti2 07 [23(]. Here we 
present a comparison of our theoretical results with data 
for Dy2Ti207. The measured magnetisation flattens off 
abruptly at the saturation value, Mexp = 10/x_b/V3, in a 
way reminiscent of the Kasteleyn transition. For temper- 
atures well below the scale set by the exchange constant, 
experiments indicate that spin ice materials are out of 
equilibrium, as expected if the actual dynamics is local. 
We therefore make a comparison at a temperature com- 
parable with the exchange energy. At such temperatures 
the behaviour of the pyrochlore lattice and BL are in- 
distinguishable. We show in Fig. [5] data for Dy2Ti207 
at Texp = 1-8K (Ref. [l2|) together with our BL results 
at a fitting temperature T ~ 0.87 Texp- This tempera- 
ture is about 1.6 times the effective nearest- neighbour 
exchange constant JeffS^/kB ~ l.llK estimated for 
Dy2Ti207 [12]. The agreement between experiment and 
our theory is clearly good, indicating that the former is 
indeed close to the transition discussed here. 

In conclusion, we propose that spin ice materials pro- 
vide a new opportunity for the study of the Kasteleyn 
transition. The 3d transition pres ented here differes from 
that in 2d, predicted to occur 2J] for spin ice materials in 
a [111] field, both in its critical properties and in the be- 
havior of correlation functions. In 2d, it takes the form of 
a commensurate-incommensurate phase transition which 



can be evidenced by a divergent peak in the correlation 
function (absent in Sd), the wavevector of which is field- 
dependent [13] , as reported in recent measurements [2^ . 
We believe that the data of Refs. [13] and [2^ provide ev- 
idence for rounded 3d Kasteleyn transitions in Dy2Ti207 
and Ho2Ti207 respectively, which could be more closely 
approached in future experiments. 
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